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Abstract 

We present a new estimate of the exponent governing the small-x behavior of the 
nonsinglet structure function g^ -n derived under the assumption that the Bjorken 
Sum Rule is valid. We use the world wide average of a s and the NNNLO QCD 
corrections to the Bjorken Sum Rule. The structure function g^ s is found to be 
clearly divergent for small x. 
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Introduction 



During the last decades many deep inelastic scattering experiments have been 
performed to elucidate the internal spin structure of the nucleon and to verify 
certain Sum Rules. One of them is the Bjorken Sum Rule [1] (BSR) which has 
been verified with an uncertainty of 8% [2] and which relates the structure 
function g\~ n to the axial vector decay constant. Perturbative corrections to 
the BSR have been calculated up to O (cvf) [3] and O (af )-corrections have 
been estimated [4] 
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with the numerical values for rif = 3. 



Above perturbative series can be used to determine a s from an experimental 
value for rf _n . In [2] it was found that such an analysis suffers from the 
poorly restrained small-x behavior of the structure function g\. Only some 
years ago it was quite common to expect g\ to be in agreement with Regge- 
theory [5] which predicts g\ ~ x x with A = . . . 0.5. Experimental data seems 
to contradict this expectation at least for the neutron structure function. A 
divergent behavior was also found by QCD analyzes of parton distributions, 
see e.g. [6,7]. There have been attempts to explain the obviously divergent 
behavior of g\ by considering perturbative effects of higher orders that may 
dominate over the leading or even next-to-leading order QCD at small x. A 
resummation of double logarithmic terms originating from ladder diagrams [8] 
led to A —0.5. It is not clear yet how large the influence of single logarithmic 
terms is and how one could perform a consistent resummation of those. 

In the first part of this letter we want to repeat the analysis from [2] and show 
how the small-x behavior prevents a conclusive extraction of a s from the BSR. 
We will furthermore show that even other uncertainties such as higher twist 
corrections or experimental errors prevent this method of determining a s from 
being competitive. 

In the second part we will make use of the fact that a s is known quite ac- 
curately from other experiments and we can therefore estimate the small-x 
exponent of g± s by assuming the BSR to be valid. 



Extracting a. 



The perturbative expansion (1) allows for a determination of a s from an ex- 
perimental value for IV This has been done in [9], where only experimental 
results from EMC, E142 and E143 have been used, and in [2]. One has to crit- 
ically analyze the way those experimental values were obtained. The problem 
is, that the structure function g x cannot be measured at one fixed Q 2 over 
the whole x-range. Therefore one has to evolve the data points to one com- 
mon scale and extrapolate into the unmeasured x-region. The first problem 
is solved using the DGLAP-equations [10] that describe the Q 2 -dependence 
of the parton distributions: 
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Table 1 

World average for higher twist terms. The last two rows show values extracted 
from experimental data for comparison, they have not been included in the average. 
Some models are able to produce Q 2 -dependent twist matrix elements. In such cases 
values for 5GeV 2 are stated. 



Method 


Year 


fP-n 
J A 




p—n 
a A,2 


ref. 


QCD sumrule 


1990 


— 0.068 ± 0.034 


0.025 ±0.012 




[15] 




1994 


-0.120 ±0.006 


0.057 ±0.003 




[16] 




1995 


-0.024 ± 0.004 


0.024 ± 0.010 




[17] 


MIT bag 


1994 


+0.035 


0.021 


0.059 


[18] 


Modified bag 


1993 




0.0052 




[19] 


Diquark 


1995 


+0.006 






[20] 


Soliton 


1998 


-0.033 


0.0003 




[21] 


Renormalon 


1996 


±0.043 






[22] 


Lattice 


2001 




0.0046 ± 0.0090 


0.034 ± 0.008 


[23] 


E142, E143, 
E154, E155 


1999 




0.003 ± 0.010 




[24] 


E143 


1996 




0.0024 ± 0.0200 


0.0311 ±0.0040 


[25] 



dAq NS (x,Q 2 ) a s (Q 2 ) NS 



d InQ 2 



, 7 i>,::®Mxs. 



AS 



d InQ 2 I AG 



a s (Q 2 
2ir 



P P 

1 qq 1 qg 

P P 

r 91 r 99 




(2) 
(3) 



with Pij being the splitting functions that have been calculated in next-to- 
leading order (NLO) [11]. The symbol ® denotes convolution with respect to 
x. With 
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p(n) _ ( e ) 



^C NS ® (±3Aq 3 + Aq 8 ) + C s ®Ai; + 2n f C G ®AG 



, (4) 



the structure functions g\ and can be related to the polarized nonsinglet 
and singlet parton distributions and gluon distributions Ag 3 8 , AS and AG, 
respectively. The corresponding Wilson coefficients have been calculated in 
NNLO [12]. The average quark charge is denned by (e 2 ) = (nj)^ 1 J2i=i e i- 

With these ingredients one can construct an algorithm that performs the 
Altarelli-Parisi evolution of the parton distributions and calculates a best 
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fit for g±(x, Q 2 ) which can then be integrated to yield the first moment 



Tt ) m=( 1 dxgf l \x,Q 2 ). (5) 

JO 

We repeated an analysis performed by the ABFR-group [6] with the NLO 
code that had been kindly provided by Stefano Forte. It parametrizes the 
input distributions as 



Aq l = r li x a *(l-xf>(l+ liX Si ) , (6) 

with i — E , G , 3 and 8. The distributions are normalized in order for i] 3 and 
i] 8 to resemble their first moments which are known from other experiments 
and can be fixed to r] 8 = 0.579 ± 0.025 [13] and r] 3 = g A /g v = 1.2670 ± 0.0035 
[14]. Of course, r/ 3 should not be fixed if one wants to test the BSR. But since 
we assume its validity and want to extract a s we may use the known value for 
gA/gv and shall not fix a s during the evolution procedure. 

We use the ABFR-procedure which fixes 5 more parameters in 4 different fit 
types, so we are left with 7 fixed and 14 free parameters (including a s ). See [6] 
for details about the fits. This procedure results in a world average for data 
from SMC, E142, E143, E154, E155 and HERMES (where only data taken at 
Q 2 > 1 GeV 2 have been used in the fit) of 



rf n = 0.1847 ±0.0035, (7) 

where systematic and statistical uncertainties from experiment have been 
added quadratically. 

This is the result of a purely perturbative NLO analysis. We may include 
higher twist corrections which are given by twist matrix elements of the form 



L 



where perturbative corrections have not been noted explicitly. The matrix 
elements a^,2 and cIa are so-called target mass corrections and arise from op- 
erators of twist 2 and 3, respectively. The higher twist corrections originate 
from an operator of twist 4 and are given by /a- Even higher twist corrections 
are suppressed with l/Q 4 and have been neglected in our analysis. There are 
many different models how to calculate higher twist corrections. An overview 
is given in Table 1. The errors produced by different models are due to some 
approximation used within the respective models and do not give a state- 
ment about the quality of the model. Therefore we did not consider them 
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Table 2 

Fit results for g\ at lowest x data points and the theoretical value for T± at the 
corresponding Q 2 . 



experiment 


X 


Q 2 /GeV 2 


gf- n) (fit) 


■ptheory 
1 1 


SMC 


0.005 


1.3 


1.1968 ±0.8516 


0.1580 ±0.0043 (match) 
±0.0162 (twist) 


SMC 


0.008 


2.1 


1.1330 ±0.6560 


0.1734 ±0.0003 (match) 
±0.0101 (twist) 


SMC 


0.014 


3.5 


1.0161 ±0.3985 


0.1806 ±0.0005 (match) 
±0.0060 (twist) 


E155 


0.015 


1.2 


0.6798 ±0.1909 


0.1548 ±0.0048 (match) 
±0.0176 (twist) 



but calculated the average value with one standard deviation as characteristic 
uncertainty. That results in a correction to the BSR at 5 GeV 2 of 



m 

w 



[a A ,2 + 4d A + 4/a) = 0.0004 ± 0.0042 . 



(9) 



The absolute size of the correction is indeed small but the uncertainty due 
to different models becomes so large that higher twist corrections cannot be 
neglected and result in a non-negligible uncertainty on a s . 

Furthermore, it has become clear that at small-x effects which are beyond 
NLO QCD become important. A resummation of double-logarithmic terms 
a™(ln 2 l/x) n [8] yields a small-x behavior for the nonsinglet structure function 
with a power-like divergence g^ s ~ x x with A ~ —0.5. This result contradicts 
Regge theory [5]. Running the evolution program again with the nonsinglet 
exponent fixed to the two extreme cases ±0.5 results in an uncertainty for 

rr n of 



p—n 



± 



0.0034 
0.0282 



(10) 



This induces a large error on a s and basically spoils the whole analysis. Our 
result for a s is summarized as follows 
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a s (M z ) = 0.1160 ± 



0.0180 
0.0043 



(small x) ± 



0.0043 
0.0054 



(twist) 



± 



0.0034 
0.0041 



(exp) 



± 



0.0031 
0.0036 



(evol) 



±0.0008 (match) ±0.0006 (g A /gv) 



The evolution error has been estimated by repeating the fit in different fit types 
and renormalization schemes and varying factorization and renormalization 
scales. The matching error results from varying the matching thresholds for 
the evolution of a s to the Z-pole from single to double quark masses, where 
we used the matching procedure from [26]. The evolution was performed with 
the three-loop /3-function. 

This result impressively states the need for more accurate experimental data 
in the small-x region and/or better theoretical predictions. Even if the small-x 
uncertainty could be significantly reduced, the experimental and higher twist 
uncertainty prevent this method of "measuring" a s from being competitive 
to other methods. However, our result is in agreement with the world wide 
average a s (M z ) = 0.1181(20) [14]. 

Also note [27] which investigates the Q 2 -dependence of certain error sources 
when extracting a s from the BSR. 



Estimating the small— re exponent 

Above result enables us to turn the question around. If a determination of a s 
suffers from the unknown small-rc behavior one can use the precise knowledge 
of a s to predict this behavior. We may ask, what behavior does g\ have to 
exhibit in order to fulfill the BSR? 

To solve this question we start from the world wide average of a s and compute 
the theoretical BSR via (1) at Q 2 given in Table 2. There, also the data points 
with the smallest x are shown, with g{~ n being obtained from our best fit 
and the experimental errors assigned to it. Note, that r^ heory is defined with 
already included higher twist corrections that were calculated as in (9). 

We now assume the difference between r^ heory and the integral over the mea- 



sured region (where we assume the extrapolation to large x to be in agreement 
with the NLO QCD fit) 




(11) 
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to be given by the power-like behavior of g\ ~ x x alone. Consequently, we 
obtain the small-x exponent by solving 

r! hcory (Q 2 ) -F^xo... 1, Q 2 ) = const r dxx x (12) 

J o 

for A, where the constant is determined from a fit through the lowest-x data 
point. The most accurate result with this method is obtained if we choose the 
data point x = 0.014, Q 2 = 3.5 GeV 2 from SMC because this point combines 
the need for a low x-value with the desire for a large Q 2 to minimize the 
influence of higher twist corrections or experimental uncertainties. We find 

A = -0.40 ±0.24 (exp) ± J;J^ (twist) ± f 5 (a s ) . 

Other uncertainties, e.g. matching, errors on qa/ov or the quark masses, 
turned out to be comparably negligible, i.e. less than 3%. 

At other data points A cannot be determined with sufficient precision, mostly 
due to higher twist corrections which become more important at such low Q 2 . 
The errors associated with the other data-points are too large to identify any 
Q 2 -dependence of the exponent A. 

An important advantage of this method is the use of the NLO code only above 
x Q , where Q 2 is relatively large. 

To check on this result we also determine the nonsinglet exponent that best 
fulfills the BSR with the help of the evolution code. We fix r/ 3 = 1.2670(35) 
and a s = 0.1181(20) in the evolution program where we simultaneously fixed 
the parameters governing the large-x behavior — being 7», 5, — to the best 
fit output values . This produces a nonsinglet exponent 

a 3 = -0.563 ±0.016, (13) 

where the given uncertainty results from experiments only. This methods ne- 
glects higher twist corrections. We may estimate the influence of the evolution 
code by varying the large-x parameters or the constants a s and i] 3 in the input 
values. This leads to only minor changes and we estimate the total uncertainty 
to be 

a 3 = -0.56 ±0.04. (14) 

The exponent a 3 agrees in LO with the nonsinglet exponent A and confirms 
our first result. As stated before, the data does not allow for a Q 2 -dependent 
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determination of the exponent, that's why it remains pointless to check for 
the agreement of A and a 3 in higher order QCD. 



A similar (only code based) method has been used in [7] with a different NLO 
evolution code. They found g 1 to be even more divergent with an error of the 
same size. 



Conclusion 

We showed that the Bjorken Sum Rule is not suited to derive a value for a s 
due to the uncertainties associated with the small-x behavior. We found 



with the major uncertainties added in quadrature. This shows the need for 
more experimental data in the small-x region. Unfortunately, the new HERMES- 
results were obtained at too low Q 2 to considerably improve this situation. 
Even if the small-x uncertainty could be significantly reduced, higher twist 
corrections are a non-negligible obstacle, they are not well enough constrained 
due to the variety of existing models. 

We provide a new QCD based estimate for the exponent governing the small-a; 
behavior of gf~ n by assuming the validity of the BSR. We find 



with errors from experiment, higher twist corrections and a s combined. This 
result was derived under the assumption that g\ is dominated by a power- 
like behavior for x < 0.014 and is confirmed by the evolution code which 
neglects higher twist corrections and agrees with A in LO only. Our result 
is in agreement with the behavior predicted by the resummation of double- 
logs and contradicts Regge-theory. The assumption of the validity of the BSR 
shrinks the width of possible exponents considerably. 

The critical point of our analysis is that we depend on the experimental value 
of gi at one data point with low x and preferably high Q 2 . Since there is no 
neutron data at our preferred data point available we calculated g[~ n from the 
fit. Repeating the analysis with different fit results does not change A beyond 
the quoted error. One could also extract from deuteron measurements but 
we considered this method to be too unprecise. 

One might also question how justified the use of a NLO code in combination 
with the NNNLO corrections to the BSR is. It has been argued in [6] that 



a s (M z ) = 0.1160 ±0.0193, 



(15) 



A = -0.40 ±0.29, 



(16) 
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this procedure is justified because of the relatively small Q 2 -span of experi- 
ments which leads to an only minor higher order uncertainty in the evolution 
procedure. Furthermore, we did rely on the NLO code only above Xq when 
determining A, in this regime one does not expect higher order QCD terms 
to be of importance due to the higher Q 2 that comes naturally with larger x. 
Besides, logarithmic terms in x do not dominate in this region. 

Of course, using the BSR in NNLO only when determining a s may change 
the absolute values in (15) but not the dominance of the small-x uncertainty, 
which is the major pillar in the second part of our analysis. However, we do not 
believe the influence of higher order uncertainties to be of major importance 
for the determination of A because of the reliability of the NLO code above 
xq. Repeating the calculation in consistency with the NLO code, i.e. with the 
NNLO BSR and the two-loop /3-function changes A to -0.52, which lies within 
the given error. 

Our analysis would profit from an improvement of estimates for higher twist 
corrections or a NNLO evolution procedure that could at least restrain the 
large-x behavior better. One might even consider to perform a two-parametric 
fit that includes higher twist corrections at the parton distribution level and 
fits higher twist corrections and the nonsinglet exponent simultaneously, al- 
though it is not a priori clear if current experimental data are already suitable 
for such an analysis. 
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